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ABSTRACT 

Context. The ESO Public Survey "VISTA Variables in the Vfa Lactea" (VVV) provides deep multi-epoch infrared observations for 
an unprecedented 562 sq. degrees of the Galactic bulge, and adjacent regions of the disk. In this survey nearly 150 new open clusters 
and cluster candidates have been discovered. 

Aims. This is the second in a series of papers about young, massive open clusters observed using the VVV survey. We present the first 
study of six recently discovered clusters. These clusters contain at least one newly discovered Wolf-Rayet (WR) star. 
Methods. Following the methodology presented in the first paper of the series, wide-field, deep JHK^ VVV observations, combined 
with new infrared spectroscopy, are employed to constrain fundamental parameters for a subset of clusters. 

Results. We affirm that the six studied stellar groups are real young (2-7 Myr) and massive (between 0.8 and 2.2 lO"* Mq) clusters. 
They are highly obscured (Ay ~ 5 - 24 mag) and compact (l-2pc). In addition to WR stars, two of the six clusters also contain at 
least one red supergiant star. We claim the discovery of 8 new WR stars, and 3 stars showing WR-like emission lines which could 
be classified WR or Olf. Preliminary analysis provides initial masses of ~30-50 Mq for the WR stars. Finally, we discuss the spiral 
structure of the Galaxy using as tracers the six new clusters together with the previously studied VVV clusters. 

Key words. Galaxy: open clusters and associations: general - open clusters and associations: individual: VVV CL009, VVV CLOU, 
VVV CL036, VVV CL073, VVV CL074, VVV CL099 - stars: massive - stars: Wolf-Rayet - infrared: stars - surveys 



* Based on observations with ISAAC, VLT, ESO (programme 087.D- 
0341A), New Technology Telescope at ESO's La Silla Observatory 
(programme 087.D-0490A) and with the Clay telescope at the Las 
Campanas Observatory (programme CN201 lA-086). Also based on 
data from the VVV survey (programme 172.B-2002). 



1. Introduction 

Borissova et al. (1201 11 1 discovered 96 new open clusters and stel- 
lar groups in the Galactic disk area that was covered by the in- 
frared (IR) VI STA V ariables in t he Vfa Lactea (VVV) survey 
(Minniti et al. 120101 Saito et al. '2012), one of six Europeam 
Southern Observatory (ESO) Public Surveys carried out with the 
new 4 m Visible and Infrared Survey Telescope for Astronomy 
(VISTA). Most of the new cluster candidates are faint and com- 
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pact (characterized by small angular sizes), highly reddened, 
and younger than 5 Myr, thus demonstrating the potential of the 
VVV survey to find new open clusters, especially those hiding 
in dusty regions. 

In the year following this discovery, our group has initiated a 
near-IR spectroscopic follow-up campaign to confirm their clus- 
ter nature, and to derive the spectral types and distances of the 
brighter cluster members. A total of 27 new VVV cluster candi- 
dates have thus far been observed, as well as ~40 other known 
cluster candidates present within the VVV area. Our final ob- 
jective is the characterization of a large sample of star clusters 
using homogeneous data and analysis methods. Our results for 
four known clusters, Danks 1, D anks 2 , RCW79 and DBS 132, 
were presented by Chene et al. (120121 hereafter Paper I). Here, 
we use the same method as in Paper I, but applied to six of the 
new VVV cluster candidates that have not been studied before 
and in which spectroscopic observations have revealed the pres- 
ence of Wolf-Rayet (WR) stars. Studies of the other clusters will 
be presented in a series of papers in preparation. 

The study of young massive clusters containing WR stars 
is a key to understand the evolution of massive stars. Typically, 
WR stars have masses of 10-25 Mq and they are believed to 
have descended from O-type stars. They exhibit strong, broad 
emission lines of He and N (WN spectral type) or He, C and O 
(WC stars). They spend -10% of their ~5Myr overall lifetime 
as WR stars (Meynet & Maeder 2005). At solar metallicity and 
in the absence of rotation, O stars with an initial mass of at least 
~25 Mq can evolve into WR stars, but most have initial masses of 
40-50 Mq, allowing them to reach the WC sequence (Crowther 
p007). Consequently, if we assume that the mass of the most 
massive star formed in a cluster is constrained by the mass of 
the cluster itself (as proposed by Weidner et al. 1201 Oi l, initial 
total cluster masses must be ~5OOM0, or even ^IOOOMq if they 
contain at least one WN or WC star, respectively (based on the 
analytical model of Weidner & Kroupa 2004). On the other hand, 
if there is no real relation between the mass of the most massive 
member and the total mass of a cluster (Eldridge I2012I I. lower 
initial total cluster masses can be accommodated. 

If main-sequence (MS) and evolved O stars are detected to- 
gether with WR stars, tentative evolutionary sequences can be 
defined. If the overall population was formed in a burst of star 
formation, WR stars have progenitors with spectral types ear- 
lier than the O stars still present in the cluster Consequently, 
the nature of WR progenitors can be constrained (e.g. Martins et 
al. |2008 in the Arches cluster) and direct connections between 
MS O stars and evolved WR stars are established. In addition, 
the identification of MS O stars is a crucial piece of informa- 
tion: it shows the position of the turn-off" of the entire population, 
providing constraints on the age (e.g. Martins et al. 2007). This 
provides upper limits on the main sequence lifetime of the pro- 
genitors of WR stars. If red supergiants are present at the same 
time as WR stars, as is the case in only two clusters (the cen- 
tral cluster and Westerlund 1) this strongly constrains the link 
between hot and cool evolved massive stars in the upper HR 
diagram, potentially revealing evolutionary sequences between 
certain WR stars and red supergiants. Obtaining information of 
clusters of different ages is necessary to establish the evolution- 
ary sequences followed by stars of different masses. Very young 
clusters have only very massive stars off" the main sequence, 
while older (i.e. 3 to 10 Myr) contain evolved O stars of lower 
and lower masses as the age increases. In short, young massive 
clusters are a gold mine to understand stellar evolution. 

Finally, once the distances to these clusters have been ob- 
tained, they can be used to analyse the Galaxy's spiral structure. 



Indeed, the unprecedented deepness and resolution of the data 
from the VVV survey allow us to extend the current sample of 
known open clusters towards the Milky Ways highly obscured 
region. The study of the VVV clusters hence opens a window 
towards unexplored parts of the Galaxy. 

We present the observational data in Section|2] summarize 
our method in Section[3] and present the details of the clusters 
in Section|4] We then discuss the likelihood that the six stellar 
groups are real clusters, as well as the Galaxy's spiral structure 
as traced by the open clusters observed in the VVV data and the 
properties of 1 1 newly discovered WR stars in Section|5]before 
we summarize this work in Section|6] 

2. Observations 

2.1. Near-IR photometry 

2.1 .1 . VVV data and photometry 

We downloaded the stacked images of the individual 2048x2048 
pixel exposures containing the six clusters discussed in this 
paper from the VISTA Science Archive (VSA) websiteQ The 
names and coordinates of our sample clusters, as well as the rel- 
evant VVV fields (see Saito et al. 120121 for more details about 
VVV data), are listed in Tableffl 

These fields were acquired in the JHKs filters in March and 
April 2010 with the VIRCAM camera mounted on the VISTA 
4 m telescope at Paranal Observatory (Emerson & Sutherland 
2010) and reduced at C AS lQ using the VIRCAM pipeHne vl.l 
(Irwin et al. |2004). The total exposure time of each of the im- 
ages was 40 s, with two images per filter on average. For a de- 
tailed description of the observing strategy, see Minniti et al. 
(120101) . During the observations the weather conditions met all 
survey constraints for seeing (< 0.8"), airmass (< 1.5), sky 
transparency (clear) and Moon distance (Minniti et al. 12010)1, 
and the quality of the data was satisfactory (i.e., not aff'ected by 
any artificial features related to problems caused by telescope 
or detector issues). JHKs false-colour images of the six clusters 
are shown in Figure [1] Additional 8 s /Tj-band images were also 
obtained to discover and monitor variable stars in these fields. 
Unfortunately, to date fewer than 10 epochs have been observed 
for our fields of interest, which is insufficient for a thorough in- 
vestigation of the presence and properties of binaries and pulsat- 
ing stars. 

Similarly to Paper I, stellar photometry was performed by 
employing the VVV-SkZ_pipeline's (Mauro et al. 2012) auto- 
mated software based on ALLFRAME (Stetson 1994), opti- 
mized for VISTA point-spread-function photometry. Extensive 
details about the many steps required to obtain photometric 
measurements are described in section 2.2 of Moni-Bidin et al. 
(1201 II) . The photometric accuracy reaches cr =0.05 mag for 
the brightest (J = 9.0 mag and/or = 9.0 mag) and the 
faintest (J = 20.0 mag and/or /Ts - 18.0 mag) stars, and reaches 
cr =0.01 mag for stars with intermediate brightnesses. 2MASS 
photometry was used for absolute ffux calibration in the J, H 
and bands using stars with 12.5 < J < 14.5 mag, 11.5 < 
H < 13 mag and 11 < < 12.5 mag. We opted to work in the 
2MASS system rather than sticking to the natural VISTA sys- 
tem, since most of the models we use in this study, as well as 
most other studies of massive clusters, employ the 2MASS pho- 
tometric system. For stars brighter than 7^5=9.0 mag we simply 
adopted the 2MASS magnitudes. 

' http://horus.roe.ac.uk/vsa/ 
" http://casu.ast.cam.ac.uk/ 
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Fig. 1. JHKs false-colour images of the six clusters. Stars marked with red circles were observed using near-IR spectrographs. 
Yellow dashed circles indicate the angular sizes of the clusters (see Section lTTT l. Coordinates are given in the J2000 system. 



2.1.2. Statistical decontamination 

Statistical field-star decontamination was performed as de- 
scribed in Paper I, using the algorithm of Bonatto & Bica (2010 ). 
This algorithm divides the full range of magnitudes and colours 



into a grid of cells with axes along Ks,(H- K^) and ( J - /rj. For 
each cell, it estimates the expected number density of member 
stars in the cluster by subtracting the respective field-star num- 
ber density. As a photometric quality constraint, the algorithm 
rejects stars with uncertainties in Kf^ and colours in excess of 
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Table 1. Sample clusters. 



VVV 


R.A. 


Dec 


I 


b 


size 


frames 


VVV field 


CL 


(J2000) 


(J2000) 


(deg) 


(deg) 


(arcmin) 






009 


11:56:03.18 


-63:18:42.19 


296.7605 


-1.1055 


15.5 X 11.6 


2xJ,2xH,12xK, 


d002 


Oil 


12:12:41.22 


-62:42:33.52 


298.5069 


-0.1710 


15.5 X 11.6 


2xJ,2xH,12xK, 


d041 


036 


14:09:03.95 


-61:15:58.14 


312.1238 


0.2125 


23.2 X 11.6 


3xJ,3xH,l8xK, 


d088 


073 


16:30:24.20 


-48:13:00.44 


335.8939 


0.1326 


23.2 X 15.5 


4xJ,4xH,40xK, 


dl05 


074 


16:32:05.86 


-47:49:30.50 


336.3733 


0.1944 


15.5 X 11.6 


2xJ,2xH,22xK, 


dl05 


099 


17:14:25.78 


-38:09:58.28 


348.7266 


0.3233 


15.5 X 11.6 


2xJ,2xH, 6xK, 


dll4 



Notes. Columns include the VVVCL name, coordinates (R.A. and Dec; J2000), the fields of view of the regions of interest, the number of frames 
in the different filters, and the name of the VVV field where the clusters are located. 



0.2 mag. The comparison fields we used for the decontamina- 
tion depended on the projected distribution of individual stars, 
clusters or clumpy extinction due to dark clouds in an image. 
Examples include a ring around the cluster or some other re- 
gion selected in its immediate vicinity. To avoid biases intro- 
duced by the subjective choice of comparison field, we com- 
pared many different geometries, positions and sizes of fields. 
Since the VVV catalogue is well populated, the differences in 
the final colour-magnitude diagram (CMD) were marginal. 

2.1.3. Completeness 

To calculate the completeness of our catalogue, we created lu- 
minosity functions for the clusters in bins of 0.5 A's-band mag- 
nitudes. For each magnitude bin, we added artificial stars to 
the cluster image within certain magnitude ranges. The relevant 
completeness fraction was calculated by recording the recovered 
fraction of artificial stars per unit input magnitude. Stellar detec- 
tion was done using DAOphot, following a similar approach as 
in the VVV-SkZ_pipeline. 

2.2. Near-IR spectroscopy 
2.2.1. Observations 

We collected spectra of most of the brightest members in 
the six clusters using the IR spectrographs and imaging cam- 
eras Sofl in long-slit mode on the New Technology Telescope 
(NTT) at ESO's La Silla Observatory, OSIRIS at the Southern 
Observatory for Astrophysical Research (SOAR) and ISAAC at 
the Very Large Telescope at Paranal (ESO), all based in Chile. 
These stars are marked with red circles in Figure [1] The instru- 
ment set-ups are summarized in Table|2] Total exposure times 
were typically 200-400 s for the brightest stars and 1200 s for 
the faintest. 

For optimal subtraction of the atmospheric OH emission 
lines, we used nodding along the slit in an ABBA pattern, i.e. 
the star was observed before (A) and after (B) a first nod along 
the slit, then at position B a second time, before returning to po- 
sition A for a final exposure. All stars were observed using two to 
five slit positions per cluster The slit position angles were cho- 
sen such that two or more stars were observed simultaneously. 
The average signal-to-noise ratio (S/N) per pixel ranges from 50 
to 150, and in a few cases S/N ~ 30. Bright stars of spectral type 
B8 to A2 were observed as a measure of the atmospheric ab- 
sorption and selected to share the same airmass as the targeted 
cluster stars during the middle of their observation. 



2.2.2. Reduction and extraction 

All reduction steps were executed using both custom-written 
Interactive Data Language (IDL) scripts and standard irai0 pro- 
cedures. Subsequent nodding observations were subtracted from 
one another to remove the bias level and sky emission lines. 
Flat fielding, spectrum extraction and wavelength calibration of 
all spectra were done in the usual way using iraf. Calibration 
lamp spectra (helium-argon and neon for OSIRIS, xenon-neon 
for Sofl and ISAAC) taken at the beginning of each night were 
used for wavelength calibration. The wavelength solution of 
each frame has an rm.s. uncertainty between ~0.3 and 0.5 pix- 
els, which corresponds to ~30 km s ' for all instruments (except 
for ISAAC: see below). Correction for atmospheric absorption 
was done using the iraf task telluric. Before running this task, 
we first subtracted a fitted Voigt profile from the Bry absorp- 
tion feature. Finally, all spectra were rectified using a low-order 
polynomial fit to a wavelength interval that was assumed to be 
pure continuum, i.e. where no absorption or emission lines were 
observed nor expected. 

Unfortunately, the ISAAC pipeline could not be used, owing 
to many problems with the data. First, we had to drop the cal- 
ibration lamp spectra, since unrecoverable sudden shifts along 
the spectral dispersion direction occurred during the night and a 
unique wavelength solution for the entire run could not be estab- 
lished. Instead, we used the OH sky lines for wavelength calibra- 
tion. Since these lines were quite faint most of the time, owing to 
the short exposure times, the resulting calibration is not optimal 
and the r.m.s. uncertainty reached ~0.8 pixels, i.e. ~50 km s"'. 
We also had to take extra care of the ISAAC spectra because 
of severe, time-variable residual fringing. The fringing pattern 
changes for every single frame and is different for the cluster 
stars and the standard star used for the correction of the atmo- 
spheric absorption. Hence it can not be removed by standard re- 
duction and extraction. The fringes are not constant along the 
spectra and hence could not be removed by simply subtracting 
a unique frequency from the Fourier transform of the spectra. 
On the other hand, the fringes in all spectra of all stars included 
in observations at a given slit position at a given time are very 
similar Luckily, almost all the stars were observed twice dur- 
ing different nights, since the weather conditions during the first 
observing run did not meet our requirements (so they had to be 
repeated by the queue observers at Paranal). In these cases, the 
following procedure was adopted. From the 3-5 spectra per slit 
position, the spectrum of an OB star was selected. A Voigt pro- 



^ IRAF is distributed by the National Optical Astronomy Observatories 
(NOAO), which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agreement 
with the U.S. National Science Foundation (NSF). 
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Fig. 2. Stellar spectra in the area of the six clusters. Hydrogen and helium recombination lines, as well as Civ, Nai and CO lines are 
marked with vertical dotted Unes. The wavelength range displayed depends on the spectrograph used. 

Table 2. Near-IR spectrograph log of observations and set-ups. 



Telescope 


Instrument 


Dates 


slit size 


Range 


Resolution 


Pixel scale 








(") 


(jj.m) 


R 


A pix"' 


VLT 8.2m 


ISAAC 


22-23 Apr 2011 


0.3 


1.98-2.43 


3000 


7.09 


NTT3.6in 


SofI 


02-03 May 2012 


0.6 


2.02-2.32 


2200 


4.66 


SOAR 4.1m 


OSIRIS 


06-07 Jun 2012 


0.42 


2.03-2.34 


3000 


3.63 



file was subtracted from the first observation (obs 1) of this star 
at each Hi, Hei and Hen line (when present). The residual, dom- 
inated by the fringes, was then subtracted from the spectra of 
all other stars observed at the same time. These new, "cleaned" 
spectra were subtracted from the corresponding spectra of the 
other observation (obs 2), and an average of the residuals was 



calculated. The averaged residuals were strongly dominated by 
the fringes of obs 2 and included very little contribution from 
the individual stellar spectra. They could hence be used to can- 
cel the fringing of all spectra of obs 2. Once the spectra of obs 2 
had been cleaned, they were used to isolate the fringes in the 
spectra of obs 1, which was, in turn, used to clean the spectra 
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of obs 1 . This procedure was performed iteratively twice to en- 
sure optimal removal of the fringes. Once the spectra of the two 
observations had been treated, they were averaged to obtain the 
final spectrum. Of course, this method of fringe removal is not 
perfect and may lead to residuals at the ~5% level, but this is suf- 
ficiently small so as to not prevent further analysis of the data. 

2.2.3. Spectral classification 

The spectra of the stars in the fields of the six clusters are shown 
in Figure|2] Spectral classification was done based on compari- 
son with available catalogues of /T-band spectra of objects with 
spectral types derived from optical studies (Morris et al. 119961 
Figer et al. 119971 Hanson et al.[l996 2005 ) as well as th e spec - 
tral catalogs of Martins et al. (2007), Crowther et al. (I2006I I. 
Liermann et al. (2009) and Mauerhan et al. (2011). The most 
prominent lines covered by our wavelength range are Br y (4- 
7) 2.1661 fim (Brackett series), Hei lines at 2.1127 fim (3p^P"- 
4s^S, triplet). Hen line at 2.188 fim (7-10), the blend of Cm; 
Niii and Oiii at 2. 1 15 fim, Niii /12. 103 yum, and Civ fines at 2.070 
and 2.079 yum. All these lines were compared with the template 
spectra from these papers. 

The equivalent widths (EWs) of the lines present in the spec- 
tral targets were measured on the continuum-normalized spectra 
using the iraf task splot. In general, we used Gaussian profile 
fitting assuming a linear background. For more complicated pro- 
files we used the deblending function or Voigt profile fitting. The 
EW uncertainties take into account the S/N of the spectra, the 
line's peak-to-continuum ratio (see Bik et al. 120051 ) and the error 
from telluric star subtraction (which was estimated at ~10-15% 
in the worst cases). The EWs of emission lines are negative by 
definition. Their values were used for qualitative estimation of 
the spectral types of OB stars using the calibrations given by 
Hanson et al. ( 1996 ) and of WR stars using Figer et al. (1997). 
In general, the spectra of most stars show Hi and Hei lines in ab- 
sorption, which is typical for O type and early-B type MS stars. 
Also, 10 stars exhibit CO, Mgi, Fe, Ca and All lines, character- 
istic of late-type stars, and are classified as K- or M-type stars 
based on comparison with templates from the spectral library of 
Rayner et al. (120091 ). Their luminosity class will be discussed 
in Section|4] We discovered a total of 11 stars showing WR- 
like emission lines, i.e. 9 WN and 2 WC types. On the other 
hand, it is possible that some of these stars are Olf stars, since 
there is no clear way to distinguish them in near infrared /T-band 
(Conti et al. 1995 ). These stars will be discussed in more detail 
in Section l53] The method of spectral classification, especially 
in the near-IR, is correct within two subtypes, so that we adopt 
this error for our estimates. 

When the S/N is sufficiently high, the luminosity class of OB 
stars is easily determined, since we can see if the Br y and Hei 
2.16 line are blended (in MS stars) or not (in supergiant stars). 
However, for many of our spectroscopic targets, such features 
are not clearly seen, due to noise. For that reason, we verify if 
the luminosity class determined for a given star is coherent with 
that obtained for the other members of its host cluster To do 
so, we measure the spectrophotometric distances obtained when 
luminosity classes V, III or I are assumed. In many cases, assum- 
ing that the cluster members are giants and/or supergiants yields 
a distance close to the edge of or outside the Milky way (Minniti 
et al. l20I II ). which makes this assumption unlikely. When all lu- 
minosity classes give plausible distances for the clusters, we use 
the distance to objects related to the clusters to take our final de- 
cision. The estimated spectral types and luminosity class of all 
spectral targets are listed in Table[3] 



3. The method 

Our methods to determine the fundamental parameters of the 
massive clusters, such as their angular sizes, extinction val- 
ues, distances, ages and masses, have already been described in 
Paper I. Here, we provide a brief summary and highlight what 
we have done diff'erently in this study. 

3.1. Angular size 

We first obtained the coordinates of the cluster centre based on 
our stellar catalogue. Starting from a first-guess value, we cal- 
culated the radial density profile (RDF) iteratively until it con- 
verged towards the best profile (i.e. with the highest and clearest 
difference between the background and the central density lev- 
els). The central coordinates of the clusters are given in Table[T] 
Once the converged RDP has been obtained, we evaluate the 
density of field stars in the vicinity of the cluster, which we de- 
fine as the background. Since we are dealing with lower-density 
clusters than in Paper I, some RDPs may not be very clear 
Hence, we systematically removed all stars from the catalogue 
that are obviously field stars. To do so, we localized the position 
of the cluster in the CMD using a decontaminated catalogue, 
assuming a cluster radius obtained from the RDP calculated us- 
ing the full catalogue. Next, we assumed that, for each magni- 
tude bin, all stars bluer than the bluest cluster stars must origi- 
nate in the field (see example in Figure[3]) and can be removed. 
This technique is optimal for highly reddened clusters and leads 
to higher-contrast RDPs. The RDPs of five of the six clusters 
discussed in this study were thus significantly improved, while 
the RDP of the sixth cluster benefited only mildly. All cluster 
RDPs are shown in Figure|4] and will be discussed individually 
in Section|4] 

To better guide our determination of the clusters' angular 
sizes, we fit a two-parameter King profile adapted to star counts 
(as in King 1966, but using star counts instead of surface bright- 
ness). The two free parameters are the projected central stellar 
density (So) and the core radius (r^). The uncertainties on the fit 
are given by the maximum and minimum values allowed within 
the observational error bars (given by Poissonian noise). The an- 
gular size of a cluster is defined as the radius where the RDP 
meets the level of the field stars. 

3.2. Cluster membership 

The present statistical decontamination approach is far from per- 
fect with respect to the determination of cluster membership. 
Indeed, the algorithm only considers the number of stars that 
are in excess (computed for a discrete distribution of colour and 
magnitude cells) within the cluster field with respect to the com- 
parison field. Since the decontamination algorithm does not take 
stellar parameters into account (not even the reddening calcu- 
lated based on the CCD) it arbitrarily selects which stars should 
be removed or retained. If a field star is kept at the expense of 
a real cluster member, it does not alter the diagnostic obtained 
from the CMD itself, since both stars have comparable colors 
and magnitudes. The ideal way to determine membership is to 
measure the radial velocities (RVs) and the the proper motions 
of the stars in the clusters region, similarly to what we presented 
in Paper I. 

Unfortunately, the VVV database currently contains obser- 
vations covering only two years, which is insufficient to deter- 
mine precise proper motions. It was not possible to construct 
proper motion histograms using an absolute proper motions cat- 
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Table 3. Position, photometry and spectral parameters of the spectroscopic targets. 



ID 


R.A. 


Dec 


J 


H 


K, 


Sp. Typ. 


E(.l - K,) 


E{H - K,) 




d 


member 




(J2()()0) 


(J2()0()) 










(mag) 


(mag) 


(mag) 


(kpc) 





WVCL009: 

1 11:56:00.21 ^63:19:01, 

2 11:56:02.37 -63:19:05, 

3 11:56:02.55 -63:19:05, 

4 11:56:03.03 -63:19:00, 

5 11:56:03.07 -63:18:54, 

6 11:56:03.78 -63:18:54, 

7 11:56:03.90 -63:18:59, 

8 11:56:04.38 -63:18:54, 

9 11:56:06.05 -63:18:59, 

WVCLOll: 

1 12:12:40.61 -62:42:31, 

2 12:12:41.13 -ei'Al-.'iO. 

3 12:12:41.36 -62:42:32, 

4 12:12:41.40 -62:42:29, 

VWCL036: 

1 14:09:00.49 -61:15:54, 

2 14:09:01.97 -61:16:02, 

3 14:09:02.37 -61:16:02, 

4 14:09:02.75 -61:15:58, 

5 14:09:02.83 -61:15:33, 

6 14:09:03.28 -61:16:05, 

7 14:09:03.37 -61:15:44, 

8 14:09:04.23 -61:16:06, 

9 14:09:04.30 -61:15:53, 
10 14:09:04.67 -61:16:13, 

WVCL073: 

1 16:30:21.57 -48:12:55, 

2 16:30:23.73 -48:13:04, 

3 16:30:23.92 ^8:13:24, 

4 16:30:23.98 ^8:13:05, 

5 16:30:24.43 ^8:13:29, 

WVCL074: 

1 16:32:05.24 -47:49:29, 

2 16:32:05.27 -47:49:14, 

3 16:32:05.46 -47:49:28, 

4 16:32:05.47 -47:49:33, 

5 16:32:05.49 ^7:49:29, 

6 16:32:05.55 ^7:49:33, 

7 16:32:05.67 ^7:49:30, 

8 16:32:05.86 ^7:49:30, 

9 16:32:05.93 -47:49:30, 

10 16:32:06.43 -47:49:26, 

11 16:32:06.99 ^7:49:21, 

VWCL099: 

1 17:14:24.26 -38:09:48, 

2 17:14:24.41 -38:09:50, 

3 17:14:24.50 -38:09:54, 

4 17:14:24.71 -38:09:49, 

5 17:14:25.42 -38:09:50, 

6 17:14:25.58 -38:09:55, 

7 17:14:25.66 -38:09:53, 

8 17:14:26.10 -38:09:49, 

9 17:14:26.34 -38:09:43, 

10 17:14:26.42 -38:09:49, 

11 17:14:26.45 -38:09:53, 

12 17:14:26.68 -38:09:45, 



05 


12.21 


12 


12 


12 


07 


98 


11.89 


11 


37 


11 


11 


19 


12.25 


11 


70 


11 


47 


72 


13.33 


12 


68 


12 


25 


63 


11.87 


11 


31 


11 


04 


44 


10.55 


9 


93 


9 


65 


49 


11.73 


11 


12 


10 


83 


44 


12.14 


11 


59 


11 


34 


26 


12.70 


11 


89 


11 


65 


30 


13.44 


13 


00 


12 


79 


71 


11.72 


10 


60 


10 


05 


18 


13.39 


12 


49 


12 


20 


86 


13.81 


12 


77 


12 


24 


01 


15.28 


12 


45 


11 


06 


80 


13.50 


9 


25 


6 


99 


90 


16.78 


13 


69 


11 


88 


94 


14.25 


10 


95 


9 


32 


48 


15.53 


12 


65 


11 


24 


87 


15.33 


12 


65 


11 


33 


62 


15.79 


12 


44 


10 


79 


59 


12.17 


11 


13 


10 


56 


78 


13.37 


10 


66 


9 


16 


99 


14.65 


12 


15 


10 


93 


84 


13.05 


11 


76 


11 


27 


90 


15.40 


12 


51 


11 


08 


68 


12.51 


11 


41 


10 


96 


48 


14.71 


11 


53 


9 


92 


53 


12.98 


10 


54 


9 


50 


13 


16.40 


13 


53 


12 


09 


25 


16.92 


13 


23 


10 


31 


10 


14.72 


11 


82 


10 


17 


99 


16.02 


13 


30 


11 


93 


80 


15.64 


12 


90 


11 


54 


19 


16.00 


13 


24 


11 


89 


00 


15.98 


13 


46 


12 


19 


48 


16.65 


13 


35 


11 


73 


92 


15.22 


12 


53 


11 


31 


46 


16.33 


12 


74 


11 


11 


05 


16.30 


13 


05 


11 


58 


61 


13.60 


11 


87 


11 


00 


76 


14.72 


13 


11 


12 


31 


59 


14.38 


12 


97 


12 


38 


35 


12.62 


10 


93 


9 


91 


40 


13.43 


11 


70 


10 


64 


97 


14.22 


12 


60 


11 


79 


72 


11.67 


10 


06 


9 


26 


58 


13.98 


12 


27 


11 


10 


81 


10.91 


9 


01 


8 


26 


87 


15.21 


12 


96 


12 


02 


23 


15.07 


13 


03 


12 


16 


05 


14.12 


12 


41 


11 


51 



B7-8V 0.18 

08-9V 0.96 

B1-2V 0.91 

BO-lVe 1.21 

08-9 V 1.02 
WN7/0f 

09-BOV 1.09 

09-BOV 0.95 

K1-2V 0.39 

B3/5V 0.74 

WN9/0If+ 1.60 

BO-IV 1.33 

B1-2V 1.71 

BO-IV 4.41 

M6-7I 5.56 

BO-1 5.05 

B2-3I 5.00 

BO-IV 4.45 

08- 9V 5.21 

09- BOV 1.76 
WN6 4.04 

BO-IV 3.93 

M7-8111 

WN9 / 04-6If+ 4.25 

BO-1 1.71 

WN7 4.70 

Ki-2ni 

08-91II 4.52 

WC8 6.18 

WN8 4.48 

BO-1 4.30 

07-8ni 4.31 

08- 09ni 4.23 
07-8ni 4.00 
M2-3ni 

WN7/04-6If+ 3.81 

M1-2I 

09V 2.78 

09- B0V 2.56 
B1-2V 2.11 

WN6+0 2.53 

WN6 2.62 

BO/09V 2.58 

WC8 1.98 

O9-B0V 3.07 

M4-5V 

M1-2V 

M1-2V 

M4-5V 



0.05 


0.12 


2.70 


yes 


0.32 


0.63 


5.52 


yes 


0.26 


0.59 


4.22 


yes 


0.47 


0.79 


5.53 


? 


0.34 


0.66 


5.25 


yes 








yes 


0.36 


0.71 


4.67 


yes 


0.29 


0.62 


5.39 


yes 


0.13 


0.25 


0.26 


no 


0.23 


0.48 


4.65 


yes 


0.45 


1.04 


10.29 


yes 


0.32 


0.87 


5.20 


? 


0.56 


1.11 


4.73 


yes 


1.46 


2.87 


1.92 


yes 


2.04 


3.62 




yes 


1.85 


3.29 


2.02 


no 


1.66 


3.25 


2.60 


yes 


1.45 


2.90 


1.80 


yes 
? 


1.75 


3.39 


1.46 


yes 


0.61 


1.15 


2.96 


? 


1.35 


2.63 


2.18 


yes 


1.32 


2.56 


2.29 


yes 








no 


1.33 


2.76 


7.47 


yes 


0.49 


1.11 


3.60 


? 


1.53 


3.06 


5.46 


yes 








no 


1.54 


2.94 


5.86 


yes 


2.53 


4.02 


1.54 


yes 


1.58 


2.92 


2.68 


yes 


1.47 


2.80 


5.45 


yes 


1.46 


2.81 


4.84 


yes 


141 


2.75 


5.07 


yes 


1.37 


2.60 


7.16 


yes 








no 


1.13 


2.48 


13.57 


yes 
? 








yes 


0.94 


1.81 


3.13 


yes 


0.83 


1.66 


5.80 


yes 


0.62 


1.37 


3.63 


no 


0.86 


1.65 


4.84 


yes 


0.91 


1.70 


6.59 


yes 


0.85 


1.68 


4.53 


yes 


0.42 


1.29 


3.34 


yes 


1.24 


2.00 


2.99 


yes 



no 
no 
no 
no 



Notes. Columns include the name of the star, its position (R.A. and Dec), /, H and A^s photometry (in mag), spectral type, E{J - K) and 
E{H - K) colour excesses (in mag), extinction (Ak) and distance (d), respectively. The final column provides our diagnostic with regards to cluster 
membership (or otherwise). 
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1 , r 




2 4 6 8 

J - K3 

Fig. 3. {J - K,) colour as a function of K, for VVV CL074 based 
on the complete photometric catalogue (black points). Only stars 
within a radius of 0.7' around the central position of the clus- 
ter are plotted. This radius was derived from the RDP calcu- 
lated from the complete catalogue. The red dots represent the 
same CMD, but based only on the decontaminated catalogue. 
The stars in the gray area are undoubtedly not part of the cluster 
and can be removed from the catalogue without any loss of in- 
formation before recalculating the RDP. Once a new and refined 
cluster radius has been determined from the RDP obtained from 
the filtered catalogue, a final statistical decontamination step was 
performed (see text). 



alogue (as done in Paper I), since the large extinction, small size 
and, in some cases, the compactness of the clusters precluded 
them from having been appropriately observed in optical low- 
resolution surveys. Indeed, based on a search of the PPMXL 
Catalog of Positions and Proper Motions on the ICRS (Roeser 
et al. 1201 Oi l, only three stars could be matched in VVV CL009, 
two in VVV CLOU, one in VVVCL036, and none in each of 
VVV CL073, VVV CL074 and VVV CL099. 

RVs were obtained for our spectroscopic targets. 
Unfortunately, the low resolution of the spectra, combined 
with the small wavelength coverage, does not allow a very 
accurate measurement of the RVs. Moreover, considering that 
the great majority of the spectra suffer from a poor S/N and/or 
fringing and/or artificial shifts in the wavelength solutions, it 
was impossible to obtain any usable measure for more than 
one or two stars per cluster. This is insufficient to provide any 
diagnostic as to a cluster's RV or to verify the membership 
probability of any of the observed stars. 

Instead, we checked the loci of the stars in the (J - /Tj) 
vs. CMD (see Figure|5]l. Of course, the results obtained this 



way can be affected by differential reddening. We attempted to 
reduce the effects of differential reddening by employing the 
reddening-free parameter 2nir = (J - H) - a(H - K^), where 
the extinction power law a - 2.14, as determined by Stead & 
Hoare ( 120091 see also Catelan et al. 1201 li for a list of several 
other reddening-free indices in the ZYJHK^ VISTA system). We 
chose this parameter to avoid the intrinsic degeneracy between 
reddening and spectral type (and since we expect to find OB 
stars in the clusters of our sample). Originally, Negueruela et al. 
(1201 W defined 0.0 < Qnir < 0.1 mag as a separating value for 
early-type stars. Subsequently, Ramirez Alegria et al. (I2012I I and 
Borissova et al. (2012) defined the separation value for OB stars 
as 0.3 < 2nir ^ 0.4 mag. Here, we prefer to plot MS isochrones 
(using red dashed lines) for an age of 5 Myr (i.e. close to the 
average of all cluster ages and young enough to clearly show 
the MS up to the highest masses). Interestingly, the choice of 
a -2.14 puts most of the MS stars close to the value Qnir = 0. 

The reddening-free CMDs (see Figure|5j can also be used 
to separate MS from pre-MS (PMS) stars (gray points), infor- 
mation which was used for our isochrone fitting and mass cal- 
culation (see below). Special care has been taken in relation to 
the cool stars in the clusters, since they could be RSGs, fore- 
ground red-clump giants or late-type MS stars. In addition to 
their position in the 2nir diagram, we plot the EW of the CO 
line, measured between 2.294 and 2.304 yum, as a function of 
spectral type. If the EW(CO) is not consistent with the relations 
of Negueruela et al. (1201 Oi l for giant and supergiant stars, we do 
not consider the star a cluster member 

We remark that, since we still do not provide an unambigu- 
ous membership diagnostic, we can not conclude wether a given 
star (projected towards a cluster) is a real member of its cluster 
or not. Here, we define a uniform approach to select the likely 
member stars used to calculate the clusters parameters. 

3.3. Extinction, spectropliotometric distance and reddening 

Our spectroscopic classification and IR observations permit di- 
rect evaluation of the clusters' extinction, distance and applica- 
ble reddening laws. Intrinsic JHK^ colours and absolute magni- 
tudes were adopted from Schmidt-Kaler (1982) for B stars and 
Martins et al. (2007) for O stars. Dust extinction at was evalu- 
ated using A^r^ = Q.6xE{J-Ks), and the distance moduli follow: 
fiQ = Ks - Mk^ - Ak^ . The standard deviation and standard er- 
ror were computed, and cases where the former is large may be 
indicative of differential reddening. 

Decontaminated ( J-H) vs. {H-K^) colour-colour diagrams 
(CCD) are presented in Figure |5] O and B stars are marked with 
blue and red points, respectively. The intrinsic colours of MS 
and giant-branch stars (Koornneef ll983l l are plotted using blue 
dashed lines. The reddening vector, plotted in orange, was de- 
termined from the calculated spectroscopic reddening. Colour- 
excess ratios, e.g., E( J - H)/E(H - K^), were tabulated for each 
cluster and to within the uncertainties the results generally match 
recent estimates of the reddening law as inferred from numer- 
ous sight-lines towards the inner Galaxy (Straizys & Lazauskaite 
2008 Majaess et al. l201 II [20121 1. except for one cluster. Sources 
located to the right and below the reddening line may exhibit ex- 
cess emission in the near-IR (IR-excess sources) and/or be PMS 
stars. 

Spectrophotometric distances to the stars for which we ob- 
tained spectra were calculated: see Table[3] For each cluster, we 
selected the stellar members for which we are confident of their 
spectral types. The resulting error is obtained using propagations 
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of uncertainties of all the values involved in the calculation of the 
distance. 

3.4. Ages 

The age can be estimated by fitting the observed CMD with non- 
rotating Geneva isochrones (Lejeune & Schaerer I20011 i. com- 
bined with PMS isochrones (Seiss et al. I2000I I. The Geneva 
isochrones are the best publicly available models for stars with 
masses greater than 10 Mq, but they are only available for solar 
metallicity. We have no direct way to determine the metallicity 
of our sample clusters, and we cannot evaluate how far from real- 
ity assuming solar values may be. On the other hand, the clusters 
are relatively close to the Sun (see Table|4]i, so solar metallicity 
is likely a reasonable assumption. 

Starting from the isochrones scaled to the previously deter- 
mined distance modulus and reddening, we apply shifts in mag- 
nitude and colour until the goodness-of-fit statistics reach a min- 
imum value (i.e. the differences in magnitude and colour of the 
observed stars from the isochrone should be minimal) for both 
the MS and PMS isochrones. For such young star clusters (i.e. < 
10 Myr) the MS isochrones are almost vertical lines in the near- 
IR and it is hard to determine the precise age, even when using 
the PMS isochrone set. For this reason, we determine an inter- 
val of possible ages based on the range of isochrones that could 
be fitted satisfactorily. The results are shown in Figure|5] where 
the PMS stars are plotted as gray points. Modelling of the spec- 
tra using CMFGEN (Hillier & Miller, 1998) is currently being done 
and will be presented in a forthcoming paper (A.-N. Chene et al., 
in prep). Once the stellar fundamental parameters of the cluster 
members have been determined precisely using these models, 
they will be used to construct HR diagrams, allowing us to com- 
pare the distribution of these stars with the latest isochrones and 
stellar tracks (including rotation) of Ekstrom et al. (120121) . 

3.5. Present day mass function and cluster masses 

We first used the cluster CMDs to select only cluster members 
on the MS. All rejected stars are either byproducts of our decon- 
tamination method (see section lZTT i or evolved (WR, RSG) stars. 
The luminosity function (obtained as described in Section lTTT l 
was converted to the present day mass function (PDMF) us- 
ing the Geneva isochrones, assuming solar metallicity and us- 
ing the intermediate ages of the clusters determined as a av- 
erage value of the age intervals as described above. For each 
cluster we established the MS turn-on point and assumed that 
all fainter objects were PMS stars. Figure |6] shows the PDMF 
of the clusters. The mass range over which the slope is calcu- 
lated is shown as vertical dashed lines: only MS objects were 
considered. An upper mass limit was set to where the PDMF 
slope starts to skew due to the effects of saturated stars, with 
the lower limit set to where the effects of incompleteness affect 
the slope. Uncertainties are also shown; these were calculated 
from the propagation of small-number stochastical uncertainties 
from Gehrels ( 1986'). Note that despite the fact that we refer to 
the cluster PDMFs and despite the young age of the clusters, 
the stellar mass function cannot be considered the true PDMF. 
Dynamical cluster evolution will result in stellar losses through 
energy equipartition, resulting in a higher average stellar mass 
and an increase in the binary fraction. This will affect the slope 
of the stellar mass function (causing it to flatten). This is sen- 
sitive to both mass and cluster lifetime (in particular for these 
clusters, which have masses of <1O^M0 and lifetimes of a few 




0.1 1 0.1 1 

R (arcmin) R (arcmin) 

Fig. 4. RDPs of the six clusters. The open circles show the RDPs 
obtained from the complete catalogue, while the black points 
are based on the filtered catalogue. The red horizontal band rep- 
resents the background level (the width is determined by the un- 
certainty on the value) and King profiles are given with uncer- 
tainties in gray. 

100 Myr; Kroupa et al. 1201 11 1 so this efl'ect must be duly con- 
sidered for these clusters. In addition, due to stellar evolution 
the number of observable massive stars with initial masses equal 
to or greater than 60 Mq decreases with time, and clusters with 
an age of 5 Myr can already have a mass function showing an 
abrupt drop at high masses. Our calculated slopes T are listed 
in TablelH with uncertainties that propagate from those on the 
individual mass bins. These slopes were used to extrapolate the 
masses of the observed main sequence stars to provide a total 
masses for the clusters. The large uncertainties on the mass are 
due to the small mass-sampling region of the PDMF slope. As a 
result the total masses that we have obtained should be treated as 
a lower limit and taken cautiously until deeper photometry can 
be analyzed. 

4. Results 

In this section, the six clusters will be discussed in more detail. 
Since all objects have one or more WR stars in their central re- 
gions, we consider them real clusters in this section, but the like- 
lihood that they could not be will be discussed in Section lSTTI 

4.1. VVVCL009 

The RDP of VVVCL009 shows clear evidence of the pres- 
ence a cluster, based on both the complete and the filtered cat- 
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alogue (open and filled circles, respectively). The central po- 
sition of the cluster associated with the optimal RDP is R.A. 
(J2000)= 11:56:03.18 and Dec (J2000) = -63:18:42.19. Our fit 
of a two-parameter King profile to the RDP reveals a back- 
ground density (Sb) level at 253 + 17 stars arcmin"^, while 
^0 = 530 + 90 stars arcmin^^ and rc = 0.17' + 0.03'. The cluster 
radius reaches 0.6' ± 0.2' (equivalent to 1.1 ± 0.5 pc using the 
distance determined below). 

Upon examination of the stars' positions in the CMD, we see 
that star 1 is located quite far ofi" towards the blue. However, it 
seems to be the effect of some differential reddening, since in 
the plot showing the reddening-free parameter, this same star is 
found along the MS. Hence, we assume that this star is a cluster 
member. As for star 4, it exhibits emission lines, but the low S/N 
of its spectrum does not enable us to determine its spectral type 
with confidence. Based on the strength of its Br y line in emis- 
sion, we identify it as a BO-1 Ve star, but we caimot distinguish 
it from a young stellar object. Hence we wiU not use it for the de- 
termination of the cluster parameters. The spectral classification 
of star 9 gives K1-2V. Its measured EW(CO) is 5.5 A, which 
corresponds to a G2-5 spectral type if the star were a giant or 
supergiant. Hence, this star is likely not a member of the cluster. 

Using stars 2, 3, 5, 7 and 8 (star 6 is omitted since uncertain- 
ties in the absolute magnitudes of WR stars are very large and 
star 1 was excluded to the determination of the A^r-value, due to 
a different reddening from the rest of the cluster members, but 
included in the determination of the distance), we obtain average 
values of Ak = 0.64 + 0.01 mag (Ay = 5.4 ± 0. 1 mag) and a dis- 
tance = 5±lkpc.TheratioEiJ-H)/EiH-Ks) = 2.11±0.05, 
which seems to be in little excess of the values for the inner 
Galaxy listed in Table 1 of Straizys & Lazauskaite (2008). This 
would be the only cluster in this study exhibiting an anomalous 
reddening law along the line of sight. This result is analogous 
to that obtained for Danks2 in Paper! Some PMS stars can 
be distinguished and are defined using Qnir = 0.30 mag and 
/Ts = 16.5 mag as diagnostic cuts. MS stars fall in the range 
-0.1 < 2nir ^ 0.1 mag. Other stars, including star 9, are found 
in 0.3 < 2nir < 0.5 mag and could be background stars that 
were not properly decontaminated, or even some contribution 
of a less reddened cluster along the same line of sight. Given the 
distance and reddening, the fit using PMS isochrones suggests an 
age near or greater than 5 Myr, while the range of ages supported 
by the MS isochrones is 4-6 Myr, which is perfectly compatible. 
The slope of VVVCL009 PDMF is F = -0.7 ± 0.2 and, using 
111 MS stars, we obtain an observed mass of 660 ± 150Mg. The 
lower limit for the total extrapolated mass is 1000 M©. 



4.2. VVV CLOU 

Based on the complete catalogue resulting from the VVV- 
SkZ_pipeline, this cluster's RDP shows a high background 
level and some over-density of stars centered on R.A. 
(J2000)= 12:12:41.22 and Dec (J2000) = -62:42:33.52. The 
profile becomes much clearer once the decontaminated cata- 
logue is used leading to a much lower background level and 
a much steeper profile. Note that the sudden decrease in the 
background level beyond a radius of ~ 3' is due to the edge 
of our photometric catalogue, hence we simply exclude this 
section of the RDP from our analysis. The parameters of the 
best-fitting King profile are 5o = 900 + 171 stars arcmin"^ and 
Tc = 0.07' + 0.01', and we also derive 5b = 66 ± 8 stars arcmin"-^. 
The cluster radius is 0.5' ± 0.1' (i.e. ~0.9pc). 



Looking at the CMD, we see that star 1 occupies a position 
that seems to be off towards the blue. On the other hand, the 
Qnir parameter places it among the MS stars, and we consider 
it as a cluster member. Star 3 was removed from the photomet- 
ric catalogue by the decontamination algorithm without leaving 
behind any other star close to its positions. This is an indication 
(but not proof) that the section of CMD space occupied by star 
3 should not contain any cluster members. Moreover, the low 
S/N of star 3's spectrum does not allow any accurate spectral 
classification and its membership can not be clearly determined. 
This leaves only stars 1 and 4 for calculation of the parameters of 
VVV CLOU, since we prefer to avoid using the WR star 2. Also, 
star 1 was excluded from the calculation of the reddening, since 
it likely suffers from local differential reddening and its inclu- 
sion would move all isochrones towards the blue despite most of 
the other points being redder. This yields A/^ ~ 1.11 ± 0.02 mag 
(i.e. Ay ~ 9.4 + 0.1 mag) and d - 5 +2kpc. We identify very few 
PMS stars; however, if we assume that the MS is located along 
Gnir = 0.0 mag, many more stars could be considered PMS 
stars. We hence are left with very poor constraints on the age of 
this cluster because of the small number of data points, and the 
cluster could either be as young as 3 Myr or as old as 7 Myr. The 
observed mass obtained using 67 MS stars is 390 + 150 Mq, and 
the lower limit for the total extrapolated mass is 950 M© when 
r= -1.0 ±0.6 is used. 

VVV CLOU is a quite small and compact cluster which 
may need to be revisited using an integral field spectrograph to 
observe more stars than can be achieved using long-slit spec- 
troscopy. Also, we note that the spectrum of VVV CLO 11-2, 
which we classified WN9, is also compatible with an Olf type. 
Later analysis will allow us to determine the real spectral type of 
this star (A.-N. Chene et al., in prep). On the other hand, we can 
already see that the distance to the star obtained is 8.0 ± 0.2 kpc 
when a OI spectral type is used, which is significantly farther 
than the distance determined for the cluster. 



4.3. VVVCL036 

This cluster's RDP based on the complete catalogue shows very 
little evidence for the presence of a cluster. Fortunately, this clus- 
ter is quite heavily reddened (see below) and the filtering of the 
field stars was done very efficiently. The cluster centre is located 
at R.A. (J2000)= 14:09:03.95 and Dec (J2000) = -61:15:58.14. 
The parameters of the best- fitting King profile are -Sq = 142 ± 
52 stars arcmin"^ and rc = 0.29' ± 0.07', and we also find 
Sb - 13 + 1 stars arcmin"^. The cluster extends beyond a ra- 
dius of 1', but the angular distribution is very patchy and our 
analysis becomes meaningless at larger radii. Hence, we adopt a 
cluster radius of around 1' (i.e. ~0.8pc). 

Stars 3 and 8 do not seem to be related to the other stars in 
the CMD. Star 3 is rather faint for its spectral type and red, and 
star 8 is very blue. However, the Qiv/s-value for star 8 places it 
among the MS stars. Therefore, we consider star 3 as a field star, 
but give no definitive diagnostic on the membership of star 8. 
As the S/N of their spectra is quite low (~30), the uncertainty 
associated with the spectral classification is larger than for our 
other spectroscopic targets, and they may be cooler than B stars. 
Star 2 is a very bright M6-71 star. Its magnitude places it 
outside the boundary of our CMD. It EW(CO) value is 33 A, 
which is compatible with its spectral classification. Considering 
this, we assume this star to be a cluster member. If confirmed 
it will be interesting case where RSG and WR stars are formed 
simultaneously in intermediate massive cluster. Star 4 is a bright 
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Fig. 5. For each cluster we present a series of three diagrams. Left: (H - /Tj) vs (J - H). The Koornneef et al. (I1983I I and Schmidt- 
Kaler ( |1982| l models for luminosity classes I and V, respectively, are plotted in blue. The O and B stars are plotted as small blue 
and red points, respectively and the reddening vector indicates the measured value of Ay- Middle: (J - K^) vs Ks- Spectroscopic 
targets are marked using red circles, but those with magnitudes or colours outside the figure boundaries are listed on the plot. The 
PMS isochrones are shown in light blue (0.2 Myr), dark blue (l.OMyr), orange (5.0 Myr) and red (lOMyr), while the two upper and 
lower limits of fitted MS isochrones are shown in light and dark green. Right: Qnir - (J - H) - 2.14(H - K^) vs A's. Spectroscopic 
targets are also marked. A 5 Myr MS isochrone is plotted using a dashed red line and the PMS stars are plotted as gray points. 



(nearly as bright as the WR star) B supergiant. The spectrum of 
star 6 has a S/N that prevents its use here. 

The OB stars are spread over a wide range of reddening 
in the CCD. This indicates that the cluster suffers from signif- 
icant differential reddening. Thus, the average value from the 
six remaining stellar members of VVV CL36, after excluding 



the WR star, is not sufficient to obtain a good determination 
of the cluster parameters. We therefore use only the CCD and 
CMD to get an average value for the reddening, and we use 
the scatter we find in our spectroscopic values as an indica- 
tion of the uncertainty in the differential reddening. This yields 
Ak - 3.0 + 0.3 mag (Ay - 24 + 3). As for the spectrophoto- 
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Fig. 5. (continued) 



metric distance, we get 2+1 kpc. Some Hii regions and YSOs 
with measured RVs are found near the cluster (Urquhart et al. 
120071 |2009] l. On average, we get a kinematic distance to the sys- 
tem of 4.3 kpc, which is in agreement with our spectroscopic dis- 
tance. PMS stars cannot be distinguished from the other stars. 
Here, we assume 2nir = 0.0 mag as the boundary, but due to 
the differential reddening the fit of the PMS isochrones does not 
yield any relevant result. On the other hand, MS stars could be 
quite well isolated using -0.25 < 2nir ^ -0.05 mag, and the 



isochrone fit yields an age range of 5-6.5 Myr The PDMF slope 
is r = -1.0 + 0,4 and leads to a lower limit for the total extrap- 
olated mass of 2200 Mq (observed mass obtained with 207 MS 
stars is 1800 + 300 Mq). 
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Fig. 5. (continued) 



4.4. VVVCL073 

The unfiltered RDP exhibits two plateaus: one that results 
in a cluster radius of around 0.5' (i.e. ~0.6pc) and another 
of around 1.0' (i.e. ~1.2pc). Considering the nature of the 
complex background (see Figure[T]), which contains obscured 
clouds, these measurements of the cluster radius may be biased. 
For this cluster, we obtain R.A. (2000) = 16:30:24.20 and Dec 
(2000) - -48 : 1 3 :00.44. The only two cluster members that could 



be observed spectroscopically are the WR stars VVVCL073-2 
and 4. The other spectroscopic targets are either near or out- 
side the cluster's radius. Moreover, stars 1 and 3 are too blue 
to be cluster members, and the EW(CO) value (18.9 A) of star 
5 would correspond to an M-type giant or supergiant if it were 
located in the cluster, while it is classified as a Kl-2 type star, 
and it is likely a field star. This leaves no OB stars from which to 
determine the fundamental cluster parameters. Fitting the CCD 
and CMD yields Ak = 2.1 mag (Ay - 19.3, with possibly sig- 
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nificant differential reddening) and an upper limit of 7 Myr for 
the cluster age. Our spectrophotometric distance of ~ 4 kpc cor- 
roborates the average kinematic distance of 4.9 kpc to the star 
forming complexes near the cluster (Russeil l2003l l. considering 
the uncertainties. VVV CL073's PDMF slope is compatible with 
the Salpeter slope. The observed mass is 800 + 200 Mq (using 52 
MS stars), and we extrapolate a lower limit for the total mass of 
1400 Mo. 

One could suppose that this apparent group of stars may be 
the result of a low-extinction window. However, we consider 
this quite unlikely, since the two WR stars are located right at 
the postulated centre, where we expect to find the most massive 
stars in a cluster Instead, we support the idea that it is either 
heavily reddened or has an unusual PDMF, with a few massive 
stars and many low-mass stars with much fainter magnitudes. 
One has to note that VVV CL073-2 could also be classified 04- 
6If. When this spectral type is assumed, the distance to the star 
(5.9 + 0.1 kpc) still corresponds to the cluster distance. 



4.5. VVVCL074 

The background of this cluster is very complex, probably be- 
cause of the presence of a dark cloud towards the North-East 
of the cluster, and the Sb level is quite uncertain. We therefore 
assume an average background value, while accounting for the 
strong variations as a function of radius. This allowed us to fit a 
King profile to the RDF, but with uncertain results. Nevertheless, 
we obtain a central position of R.A. (J2000) = 16:32:05.86 and 
Dec (J2000) = -47:49:30.50, Sr = 150 + 20 stars arcmin-^, 
So = 880+ 114 stars arcmin-2^jj(j^^ = 0.12' + 0.02'. The cluster 
radius is 0.6' + 0.2' (i.e. 1.1 + 0.5 pc). 

The MS observed in the CMD of this cluster appears quite 
broad. The CCD shows evidence of some (but not too severe) 
differential reddening. Interestingly, the reddening-free CMD 
also shows significant scatter in the Qnir distribution, almost as 
if two sequences were present. Among our spectroscopic targets, 
the OB stars are found at -0.35 < Qnir ^ -0.15 mag, while 
cool stars and star 9 (WN7 or 04-61f) attain higher values of 
Qnir- One has to note that the classification of VVV CL074-9 
is still ambiguous, since it could be WN7 or 04-6If. The dis- 
tance to that star would be 6.8 + 0.3 kpc if it is an O super- 
giant star, which is within the cluster considering the uncertain- 
ties. Two cool stars are present in the CMD (stars 8 and 11) 
which could be cluster RSGs. While this is plausible for star 
11, which has EW(CO)= 19.2A and is classified as M1-2I, it 
is not for star 8, which has EW(CO) = 14.7 A and is classified 
as M3111. Consequently, only star 1 1 is considered a likely clus- 
ter member. Using all OB stars, i.e. stars 1, 4, 5, 6 and 7, we 
get Ak = 2.78 + 0.01 mag (Ay = 23.4 + 0.1 mag) and a dis- 
tance of 6 + 1 kpc. The luminosity class 111 has been attributed 
to the OB stars, since their spectra show that they are clearly 
not MS stars, and because their distance would have been very 
close to the edge of the Milky Way if they were supergiants. 
Moreover, this spectrophotometric distance is compatible with 
a kinematic distance of 4.9 kpc (determined using the same re- 
gions as for VVV CL073). This cluster does not seem to contain 
many PMS stars, if any. However, it is possible that PMS stars 
cannot be clearly isolated in the Qnir diagram because of the 
large uncertainties described above. The range of possible ages 
for this cluster is 4-5.5 Myr The PDMF slope of VVV CL074 is 
highly unreliable, due to the small mass-sampling region, hence 
the Salpeter slope (r=-1.35) was imposed. The observed mass 




Fig. 6. PDMFs of the VVV CL009, VVV CLOU, VVV CL036, 
VVV CL073, VVV CL074, VVV CL099. The slope is calcu- 
lated between the mass limits denoted by vertical dashed lines. 
Overplotted as a red dashed line is the least-squares fitted slope. 



is 1900 ± 200 Mq (using 170 MS stars), and the lower limit of 
the total cluster mass is 2500 Mq. 



Table 4. Summary of the physical cluster parameters. 



VVV 


radius 


Av 


d 


age 


r 


Mass 


CL 


(') 


(mag) 


(kpc) 


(Myr) 




(Mo) 


009 


0.6 ±0.2 


5.4 ±0.1 


5 ± 1 


4-6 


-0.7 ±0.2 


660±150 


on 


0.5 ±0.1 


9.4 ±0.1 


5±2 


3-7 


-1.0 ±0.6 


390±15O 


036 


> 1. 


24 ±3 


2± 1 


5-7 


-1.0 ±0.4 


1800±300 


073 


0.5 ± 0.5 


~ 19.3 


~ 4 


< 7 


-1.33 ±0.5 


800±200 


074 


0.6 ± 0.2 


23.4 ±0.1 


6± 1 


4-6 


Salpeter 


1900±200 


099 


0.8 ± 0.2 


15.2 ±0.2 


4± 1 


4-6 


-1.2 ±0.2 


1000±100 



Notes. Columns include the VVV CL name, the cluster radius (r^), its 
reddening, distance, age, the slope of the PDMF ("Salpeter" means that 
the slope has been imposed) and the observed mass. 



4.6. VVVCL099 

The RDP of VVV CL099 is not like the RDP of a perfectly 
symmetric and nicely shaped cluster Nevertheless, we could 
obtain a central position of R.A. (J2000) = 17:14:25.78, Dec 
(J2000) = -38:09:58.28. Filtering of the catalogue enabled us 
to obtain a better-defined profile, but our fit of a King profile 
was not fully satisfactory. We can still use the values we ob- 
tained as first estimates, with - 400 + 30 stars arcmin"^, 
5o = 638 + 67 stars arcmin^^ and = 0.37' + 0.08'. The cluster 
radius is 0.8' + 0.2' (i.e. 1.2 + 0.5 pc). 

Stars 3, 9, 10 and 1 1 can be rejected based on their position 
in the CMD, combined with their spectral types. Moreover, the 
cool stars 9, 10 and 11 cannot be classified as giants or super- 
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giants, since their EW(CO) values (15.8, 9.4 and 11.6 A, respec- 
tively) would not correspond to their spectral types (M2-4V) in 
that case. Similarly, star 12 has EW(CO) = 16.8 A and is classi- 
fied as M2^V too. Other stars, excluding the WR stars stars 4, 
5 and 7, can be used to derive the cluster parameters. We obtain 
Ak - 1.79 + 0.02 mag {Ay - 15.2 + 0.2 mag) and a distance of 
4+1 kpc. Stars with 2nir - 0.49 mag are defined as PMS stars. 
Interestingly, while our OB spectroscopic targets are close to the 
MS isochrone in the reddening-free CMD, the rest of the MS 
seems shifted towards higher 2nir values. However, the mod- 
elled MS can be aligned to the observed sequence if star 8 is ig- 
nored for the determination of the cluster parameters. The fit of 
the PMS isochrones coincides with ages older than 5 Myr, while 
the fit of MS isochrones yields an age between 4 and 6 Myr, 
which are mutually consistent. The PDMF slope is compatible 
to the Salpater slope and the observed mass, determined using 
168 MS stars, is 1000 + 100 Mq. This gives a lower limit for the 
total extrapolated mass of 1600 M©. 

5. Discussion 

5. 7 . On the likelihood of these objects being real clusters 

The new stellar groups discovered by Borissova et al. ( 1201 11 1 
exhibit significant spatial stellar overdensities. However, only a 
follow-up study can confirm if these are real clusters or sim- 
ple asterisms. This study covers six of these groups in which 
spectroscopic observations have revealed the presence of at least 
one WR star, sometimes associated with as many as 5-9 OB 
stars. The presence of this many OB and WR stars inside such 
a small angular radius is an indication of clustering. It cannot 
be completely excluded that O stars are formed outside a cluster 
environment, since a certain fraction (~ 4%) are observed to be 
apparently isolated (de Wit et al. 2005). One could then wonder 
if a random distribution of field O stars can produce an align- 
ment which may give the impression of clustering. However, it 
seems that this same fraction of O stars are formed in clusters 
which contain no other massive stars, giving the impression of 
being isolated, when stars are sampled randomly from a standard 
PDMF and a standard cluster mass function (Parker & Goodwin 
I2007i ). Finally, field O stars ejected from nearby clusters (run- 
aways) are the only real isolated stars and should always exhibit 
peculiarly high proper motions and/or RVs (e.g. Gvaramadze et 
al. 2012). 

In four of our sample clusters, 5-9 OB stars have spectropho- 
tometric distances of up to 1 kpc, i.e. within the uncertainties 
of the distances determined for their associated clusters (note 
that these uncertainties are dominated by the uncertainties in the 
spectral classification). Moreover, all OB stars in a given cluster 
are found along the same MS in the CMDs and are never sepa- 
rated by an angular distance greater than ~ 1'. Considering the 
distances to the clusters, this means that all OB stars are confined 
within 1-2 pc. As for the other two clusters (VVV CLOU and 
VVV CL073), they contain one B and one WR star, and two WR 
stars, respectively, and are confined within less than 1 pc. All six 
clusters have WR stars near their centers, i.e. where massive stars 
are preferentially found when clusters are compact (Bate 2002; 
Allison et al. 2009). Can a random distribution of massive field 
stars mimic what seem to be typical cases of massive, compact 
clusters? 

Let us, for a moment, assume that O stars can form in iso- 
lation, and have a look at how O- and B-type stars would be 
distributed in the Galactic plane when one randomly picks stars 
from a standard PDMF along the lines of sight to our sample 



clusters. Using the Besan^on model of Galactic stellar popula- 
tion synthesis (Robin et al. '2003), we looked at the distribu- 
tion of stars with spectral types between 07 and B5 (more or 
less the spectral types observed in this study) with distances be- 
tween and 50 kpc. Because we are studying clusters at low 
Galactic latitudes, we could in principle use highly diffuse ex- 
tinction (i.e. 0.7 mag kpc"') for our model, but in case we are 
dealing with windows of lower extinction, we used normal to 
low-level diffuse extinction and look for the upper limit of the 
number of OB stars per unit surface area. The simulation was 
run using 35 patches of sky, covering 1 1 arcmin^ each, around 
each of the clusters' coordinates. On average, we obtained a to- 
tal of 0.2+0.1 O- and B-type stars arcmin"^, which is well below 
what is needed to reach the densities observed (i.e. values of 2.5- 
8 O and B stars arcmin"^, depending on the cluster). We hence 
affirm that the six stellar groups studied here are not asterisms, 
but real young, massive and compact open clusters. 

5.2. Milky Way structure 

In this section we discuss the spiral structure in the fourth quad- 
rant of the Milky Way based on the location of objects related 
to very young populations (< 10 Myr) which can be considered 
reliable tracers of spiral structure. In Figure|7]we show a hybrid 
map of this region of the Galaxy, plotting our results in com- 
bination with data from previous papers. We have included the 
embedded clusters studied in this paper, as well as those from 
Paper I and preliminary results from Baume et al. ( 1201 Oi l. We 
also include optical results from studies of young open clus- 
ters and young populations near the Carina region (Vazquez et 
al. (20051 Baume et al. 2009; Carraro & Cos ta |2009) and long- 
period Cepheids (see Majaess et al. I2009I I. We also plot the 
Hn-region compilation of Hou et al. ( 120091 1 and the four-arm 
Galactic model of Vallee ( I2008I I. We adopt 8 kpc as the distance 
to the Galactic Centre. 

Considering only the clusters and young populations, three 
clear groups can be noticed: a trace defined mainly by clusters 
detected in the optical regime, another parallel trace of clusters 
studied in the IR based on VVV data and finally a group of em- 
bedded clusters closer to the Galactic Centre. The first two traces 
follow the spiral structure indicated by Hii regions and Cepheids, 
but they seem offset with respect to each other In addition, none 
match the Vallee model, which places the Carina arm between 
them. This offset is real and cannot have been caused by sys- 
tematic effects, since a systematic error on the distance to the 
clusters (the only source of uncertainties in this analysis) would 
not preserve the spiral shape. There are two interpretations for 
the difference between the traces of the clusters detected in the 
optical and those discovered using the VVV data. First, both may 
trace two different Galactic arms (Carina and Centaurus, as in- 
dicated in Paper I) and the Vallee model would consequently 
need different parameters (pitch angle and/or interarm distance) 
to fit the fourth quadrant. Second, if we assume that the Vallee 
model is right, both traces may reveal the internal structure of 
only one arm (the Carina arm), as can be seen in other galaxies 
(e.g. M51). Clusters discovered in the VVV data would be em- 
bedded in dark clouds behind the young (optical) populations at 
the front end of the arm. Figure[8] illustrates the projected clus- 
ter distribution along the Carina arm, reinforcing the second in- 
terpretation. However, there is evidence of discrepancies among 
cluster distances obtained using different methods and this could 
affect the interpretation of the data (see Moises et al. l201 ll for a 
full discussion). 
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Fig. 7. Comparison of the observations to the Vallee (I2008I I 
model with four arms in the fourth quadrant of the Galaxy. 
We assumed a distance to the Galacic Centre of Rq - 8.0 kpc. 
Cepheid distances (Majaess et al. 2009) are plotted as green 
points, open clusters and young groups from optical iUBVI) 
data (from several studies) are plotted as blue points, while red 
points show near-lR clusters studied using VVV data (Paper 1, 
this study and Baume et al. 2010). Red circles mark the Hii trac- 
ers (Hou et al. ,2009) and the yellow dot marks the position of 
the Sun. 



5.3. The 11 newly discovered WR stars 

Here we provide a summary of the 11 stars showing WR- 
like emission lines that were discovered in our sample clus- 
ters. These stars are VVVCL009-6 (WN7/Of), VVVCLOll- 

2 (Olf or WN9), VVV CL036-9 (WN6), VVV CL073-2 and - 
4 (04-61f or WN9 and WN7, respectively), VVV CL074-2, - 

3 and -9 (WC8, WN8 and 04-61f or WN7, respectively), and 
VVV CL099-4, -5 and -7 (WN6, WN6 and WC8, respectively). 
The three stars which were classified with an ambiguous spec- 
tral type due to the similarities of the spectral lines of Of and 
WR stars in the /T-band will be studied further in an upcoming 
paper (A.-N. Chene et al., in prep), and will be consider as WR 
for the rest of this sections. We remark that if all these stars are 
O supergiants, VVV CL073-2 and VVV CL074-9 would be both 
close to and within their host cluster, respectively. Following the 
nomenclature used in the Vlll* catalogue ofWR starfl we can 
attribute a WR number to them. Table|5]includes some additional 
information that is not already included in Table[3] 

The 1 1 WR stars have late spectral type, which corresponds 
to the Galactic population (cf Crowther 2007). Also, none of 
them have a WNLh spectral type (like in Danks 2, see Davies 
et al. 120 121 and Paper I), imposing a minimal age of 2-3 Myr for 
our six sampled clusters. According to the population of O stars 
found in the clusters, we can deduce that the precursors of the 
WR stars presented here were intermediate MS O stars (i.e. typ- 
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Fig. 8. Left: Distance from the center of the Carina arm (as mod- 
eled by Valle "2008) of the clusters shown in Figure|2]as a func- 
tion of their position along the arm. The Centaurus (over) and the 
Crux (under) arms are also plotted as slopes. Right: Distribution 
of the clusters with respect to the center of the Carina arm. The 
distribution of the clusters from optical data is plotted in blue, 
and the distribution of the clusters studied using VVV is plotted 
in red. 

Table 5. Additional information about the new WR stars 



name 


WR 


/ 


b 


d 


VVV 


num. 


(deg) 


(deg) 


(kpc) 




VVVCL009 










6 


45d 


296.7623 


-1.1086 


5 + 1 


VVV CLOU 










2 


46ae 


298.5065 


-0.1702 


5±2 


VVVCL036 










9 


60aa 


312.1248 


0.2135 


2+1 


VVVCL073 










2 


75bda 


335.8921 


0.1327 


^ 4 


4 


75bdb 


335.8925 


0.1321 


tt 


VVVCL074 










2 


75bfa 


336.3755 


0.1987 


6± 1 


3 


75bfb 


336.3731 


0.1957 


It 


9 


75bfc 


336.3734 


0.1942 


It 


VVVCL099 










4 


84f 


348.7265 


0.3276 


4± 1 


5 


84g 


348.7277 


0.3256 


It 


7 


84h 


348.7274 


0.3244 


It 



ically 05-8V), corresponding to initial masses ~3O-5OM0. Of 
course, more details will be obtained after the spectral analysis 
are presented in the next study. 



6. Summary 

This paper presents the first study of six new, young, massive 
clusters from the VVV catalog (Borissova et al. 1201 II) . These 
clusters contain at least one WR (the VVV CLO 11-2 could be 
04-61f star) star, and one or more O- and B- type stars. They are 
highly reddened (Ay ~ 5 - 24 mag) and compact (~ 1 - 2 pc). 
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Their age is less than 7 Myr and they have total masses between 
0.8 and 2.2 10^ Mq. It is highly improbable that the six stellar 
groups are not real clusters but instead a simple alignment of O- 
and B- type stars along the line-of-sight. 

We report the discovery of a total of 8 WR stars, and 3 stars 
showing WR-like emission line which could be either WNL or 
Olf stars. From a preliminary study of the clusters' population 
of WR and O stars, we estimate the initial mass of the most mas- 
sive star in the sampled clusters was 30-50 Mq and determine an 
minimum age limit of 2-3 Myr for all six clusters. The clusters 
VVV CL036 and VVV CL074 contain one or more WR star and 
probably a least one M supergiant star This cohabitation will be 
discussed in an upcoming publication. 

Finally, the six sampled clusters, joined with the previously 
studied VVV clusters, were used to trace the spiral structure of 
the Milky Way. Interestingly, the VVV clusters follow the spi- 
ral structure indicated by Hii regions and Cepheid stars, but are 
clearly offset with respect to the open clusters discovered in opti- 
cal light. The two possible explanations are that both optical and 
VVV clusters trace different Galactic arms, or that they trace 
different parts of the same arm, in this case, the Carina arm. 
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